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X-ray bursts as a probe of the corona: 
the case of XRB 4U 1636-536 
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ABSTRACT 

To investigate the possible cooling of the corona by soft X-rays bursts, we have studied 
114 bursts embedded in the known X-ray evolution of 4U 1636-536. We have grouped 
these bursts according to the ratio of the flux in the 1.5-12 keV band with respect 
to that in the 15-50 keV band, as monitored by RXTE/ASM and Swift/BAT, respec- 
tively. We have detected a shortage at hard X-rays while bursting. This provides hints 
for a corona cooling process driven by soft X-rays fed by the bursts that occurred on 
the surface of neutron star. The flux shortage at 30-50 keV has a time lag of 2.4±1.5 
seconds with respect to that at 2-10 keV, which is comparable to that of 0.7±0.5 
seconds reported in bursts of IGR 17473-2721. We comment on the possible origin of 
these phenomena and on the implications for the models on the location of the corona. 

Key words: stars: corona - stars: neutron - X-rays: individual (4U 1636-536) - X-rays: 
binaries 



1 INTRODUCTION 

In low-mass X-ray binaries (LMXBs), a neutron star (NS) or 
a black hole (BH) accretes matter from a Roche-lobe-filled 
donor star and releases the gravitational potential-energy 
mostly in X-rays, either persistently, or violently in form 
of the outbursts (examples of outbursts have been shown 
in lower panel of Fig.l). The outbursts usually last from 
a few days to several months due to a sudden increase in 
accretion rate and va ry in intensity by a factor of J5 100 
iJGallowav et al.l 120081 ). During the evolution of an out- 
burst, the system usually experiences a sequence of spectral 
states, which can be simplified into a high/soft (HSS) and 
a low/hard state (LHS). For some X-ray binaries (XRBs), 
during the decay of the outburst the spectrum becomes hard 
again at a lower luminosity level, forming a so-called hys- 
teresis characterized with a lagging LHS. Such hysteresis 
was s o far observed in a NS XRB only for IGR J17473- 
2721 (jZhang et al.ll2009l V The persistent emission of the 
outbursts in NS LMXBS was described by the sum of a 
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black body component plus a thermal Comptonization spec- 
trum (jTitarchukl 1 19941 ) . The blackbody component (at sev- 
eral keV) was attributed to the emission from the NS sur- 
face and the innermost accretion disc. The Comptonized 
component (at tens to hundreds of keV) was thought to 
originate from a hot corona. However, the location and for- 
mation process of the latter is still poorly known. Accord- 
ing to the proposed corona location, models for NS XRBs 
are divided into the Eastern model, for which the coron a 
is in the vicinity of the NS surface (jMitsuda et al.lll989l ). 
and t he Western model, for which the corona surrounds the 
disc (IWhite et al. 198o). Regarding the corona f ormation, 
evaporation (Mevcr ct al. 1994; Esin ct al. 1997]; iLiu et al.l 



120071: iFrank et al.l 120021') and magnetic re-connection mod- 



els (|Zhang et al.ll2000l : IZhang||2007l : iMaver et al.ll2007l ') have 
been put forward. 

Type-I X-ray bursts are thermonuclear explosions on 
the surface of neutron stars other than releasing the gravi- 
tational potential-energy. They manifest as a rapid increase 
in X-ray intensity, many times brighter than the persis- 
tent emissions, followed by distinct spectral softening during 
bursts decay, which is the result of the cooling of the neutron 
star photospheres. Unlike the timescale of days to months 
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of the outbursts, the timescale of typical type-I bursts is 
tens to hundreds of seconds. Type-I X-ray bursts, unique to 
NS XRBs, c an be used to probe the corona. For instance, 
IChen et al.l (|2012l ) observed IGR J17473-2721 during the 
2008 outburst |j They found that for bursts embedded in 
the preceding LHS, the 2-10 keV and 30-50 keV fluxes are 
anti-correlated, and that there is an apparent time lag of 
abou t 0.7 ± 0.5 seconds present in the latter (IChen et al.l 
I2OI2D . 

Here we report on LMXB 4U 1636 -536 (also known a s 
V801 Ara), discovered with KOSO-8 (jSwank et al.lll976h . 
and subsequently st udied with SAS-3 , Hakucho, Tenma, 
EXISAT, and RXTE (JLewin et al.lll987l '). The spectral clas- 
sification wit h the color-co lor diagram (CCD) has been re - 
ported by Hasinger et all ([1989) and ISchulz et all (|l989l) . 
In the X-ray CCD or hardness-intensity diagram (HID), 
4U 1636-536 traces a U-shape or C-shape as a typical atoll 
source (detailed CCD and HID were shown in IZhang et al.l 
l|201lh ). ( The type-I bursts in Z-sources are relatively rare 
compared to in atoll sources due to large accretion rate and 
the presence of a stro ng magnetic field. (jLewin et al.lll993l : 
ICallowav et all 120081 ) ) The orbit and spin periods o f 4U 
1636-5 36 are about 3.8 hours (jvan Paradiis et al.lll99d ) and 
581Hz l|Strohmaver et al.lll998alb|), res pectively. Its distance 
was estimated as ~6 kpc (|Gallowav et a l. 2006') during the 
type-I bursts with photospheric radius expansion (PRE). 

The rest of this work is organized as follows. Section 
2 introduces the data and the observations. Section 3 de- 
scribes the findings of additional hints for an anti-correlation 
between soft and hard X-rays in bursts. Under these find- 
ings. Section 4 describes the possible constraints upon the 
disk/corona configuration during state transitions. A discus- 
sion on the implications of these results is offered in Section 
5. 
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Figure 1. Upper panel; the long-term monitoring of 4U 1636-536 
at soft X-rays (1.5-12 keV) by ASM (black), and at hard X-rays 
(15-50 keV) by BAT (blue). The red lines are the type-I bursts 
occurred during the evolution of the source. Lower panel: example 
of the outbursts recorded by ASM (black) and BAT (blue). The 
locations of the bursts are symbolized by asterisks and classified 
into groups A and B. 



2 OBSERVATIONS AND DATA ANALYSIS 

4U 1636-536 has been monitored by the All Sky Monitor 
(ASM) of Rossi X-ray Timing Explorer (RXTE) at 1.5-12 
keV since 1996. Here, we consider the RXTE observations in 
MJD 53452 - 55847, during which monitoring at 15-50 keV is 
also available from the Burst Alert Telescope (BAT) of Swift 
(see Fig.l). [j On top of Fig. [TJ we mark 126 type-I bursts 
embedded in the evolution of the source. Among them, 12 
PRE bursts are not included in the following analysis since 
their flat top peaks make it hard to co-align the individual 
burst according to the flux evolution. For the remaining sam- 
ple of 114 burst events, we extracted their li ghtcurve from 
the P roportional Counter Array (PCA) data l|jahoda et al.l 
Il996l ). We utilize the best-calibrated PCU2 only. The anal- 
ysis of the PCA data is performed using HEAsoft V. 6.12. 
Light curves in bins of 1 second are extracted from the sci- 
ence event data of E_125u_64M_l_s. The persistent emission 



1 The outburst in IGR J17473-2721 had a long-lived preceding 
low/hard state (lasting about two months), this stable plateau 
represents a significant difference in comparison to the sharp peak 
usually seen in atoll sources, like Aql X-l(Zhang et al, 2009). 
^ For more information about ASM and BAT, see; 
[http;//heas arc. gsfc.nasa.gov/xte_ weather/^ 
and http;//swift. gsfc.nasa.gov/docs/swift /results/transients/ 



is taken at a window of 30 s prior to each burst, and sub- 
tracted off to have the net burst lightcurve. The dead time 
correction is made under the canonical criterion described 
at the HEASARC Web sit^, with the standard lb data as 
a reference. Lightcurves for the bursts are co-aligned and 
combined according to the peak flux of each burst in the 
2-10 keV band. 



3 RESULTS 

The 4U 1636-536 outbursts adopted in our analysis are 
shown in Fig. 1. The average peak fluxes of the outbursts 
are derived as 12 ct/s at 1.2-12 keV from the total ASM 
lightcurve, and 0.015 ct/cm^/s at 15-50 keV from the total 
BAT lightcurve. 

A zoom-in view of a few outbursts is shown in the 
bottom panel of Fig. 1. The spectrum evolves from the 
low/hard state to the high /soft state during an outburst 
(|Lewin et al.ll 19931 : ICallowav et al.ll2008l '). Therefore, the al- 
ternately hard/soft state can be regarded as numerous sepa- 
rated outbursts, in which the hard X-ray peaks lead the soft 
X-ray ones. We subdivided the ASM and BAT lightcurves 



^ Please see the website: http://heasarc.nasa.gov/docs/xte/recipes 
/pca_deadtime.html 
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Figure 2. The time delay distribution for the outbursts of 4U 
1636-536. The time delay refers to the time in which the soft X- 
ray peak (monitored by ASM) lags the hard X-ray one (monitored 
by BAT). 
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Figure 3. Lightcurves derived from PCA for bursts combined 
in group A: 1 s-bin lightcurve at 2—10 keV (top panel), 8 s-bin 
lightcurve at 30—50 keV (bottom panel). 



into 30-days segments, covering in each both the individual 
hard and lagging soft states, and calculated the cross cor- 
relation between them |j and then utilized a throwing-dot 
method, further described below, to determine the time lag 
and its error. The distribution of the time lag for the out- 
bursts is shown in Fig. 2, presenting an average of about 12.7 
±0.2 days. This time lag is far beyond zero which means 
the hard/soft states can be distinguished clearly, and hence 
the type-I bursts can be well classified accordingly. 

We subdivided the bursts into two groups. Group A is 
composed of 78 bursts, which are located around the out- 
burst peaks in the ASM lightcurve and the each valley of 
the BAT lightcurve. Group B encloses 36 bursts embedded 
in the preceding LHS, with a relatively high persistent hard 
X-ray flux, and hence occurring in hard states of the persis- 
tent emission. The average hardness ratios (BAT flux/ ASM 
flux) turn out to be (0.4 ±0.1) x 10"^ and (1.4 ± 0.4) xlO"^ 
for groups A and B, respectively. Such a classification allows 
investigating the possibly different outburst evolution. 

We combined the 2-10 keV and the 30-50 keV 
lightcurves for group A and averaged them in each time bin, 
taking the 2-10 keV peak time of each burst as reference. 

^ The cross correlation was calculated with the software 
CROSSCRO, and detailed information is available at 
|http://heasarc.gsfc.nasa.gov/ftools/fhelp/crosscor.txt| 
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Figure 4. Lightcurves derived from PCA for bursts combined 
in group B: 1 s-bin lightcurve at 2-10 keV (top panel), 8 s-bin 
lightcurve (middle panel), and revised lightcurve (bottom panel) 
at 30-50 keV. 
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Figure 5. Evolution for the average burst fluxes in the 2-10 keV 
band against the ones at the 30-50 keV band. The 8 s-bin average 
burst fluxes are derived from combining the bursts enclosed in 
group B, which corresponds to the up and bottom panel in Fig. 4. 
The line shows a linear fit to the data. 



The resulting averaged lightcurves are produced in bins of 
1 s at 2-10 keV and 8 s at 30-50 keV. The results are shown 
in Fig. 3. We notice that the 30-50 keV flux close to zero 
of the time axis shows a weak (statistically non-signiflcant) 
excess co-aligned with the soft peak at 2-10 keV. If it is 
not a statistical fluctuation, this excess may come from the 
reprocessing of the burst seed photons by the corona. 

We have obtained the averaged lightcurves at 2-10 keV 
and 30-50 keV for group B as well (see the top and middle 
panels of Fig. 4) . There exists a clear shortage in the 30-50 
keV flux, a phenomenon similar to that discovered in the 
preceding LHS of the 2008 outburst of IGR J17473-2721 
(Chen et al. 2012). A constant fit to this lightcurve results 
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4 STATE TRANSITION 

During the outburst of an XRB, the transition from LHS 
to HSS is accompanied with a large suppression of the hard 
X-rays, which is usually understood as the corona cooling 
by the shower of soft X-rays from the disk. Since the type-I 
bursts cool the hot corona in a similar manner, and given 
that the luminosity and location of the type-I bursts are well 
known, we may use the shortage observed in hard X-rays 
while bursting to constrain the corona height with respect 
to the disk at the time of the state transition. Consider the 
following definitions of the energy densities encounted by 
the corona: 



(-^burst 



L 



47rc(r-c2-f ft2)' 



(1) 



Figure 6. The cross correlation between fluxes of 2-10 keV and 
30—50 ke V for bursts combined in group B . The fluxes are derived 
in bins of 1 s. The curve shows a Gaussian fit to the data for esti- 
mating the time lag in these two energy bands. The y-axis denotes 
the cross-correlation factor derived between the lightcurves in 2- 
10 and 30-50 kcV bands. 



in a X of 49.40 under 25 dof, suggesting a significance of 3(t 
for the shortagelj 

With the 2-10 keV and the 30-50 keV lightcurves for 
group B, we investigated the correlation along the burst evo- 
lution between the hard and soft X-rays (Fig. 4) . Lightcurves 
in both energy bands, which corresponds to the upper and 
lower panel of Fig. 4, were produced with 8 s bins; the cor- 
responding fluxes are plotted in Fig. 5. A linear fit to this 
flux-flux diagram gives a slope of — (7.3±1.6)xl0~* under 
a reduced x^ of 0.7 (15 dof), and a correlation coefficient 
of —0.8. The possible time delay between the two energy 
bands was investigated with 1 s-bin lightcurves as well|j We 
utilized a cross correlation and a throwing-dot method to es- 
timate the time lag and its error, and derived a value of 2.4 
± 1.5 seconds. In practice, the throwing-dot method is that 
for bursts in group B we combined the 1 s-bin lightcurves 
at 2-10 keV and 30-50 keV, respectively, and conservatively 
subtracted off the excess in the 30-50 keV band, although 
it has little influence on the time delay. Then we sampled 
the lightcurve by assuming that the flux for each bin has 
a Gaussian distribution, and estimated the time delay with 
a cross-correlation method. By sampling the lightcurve a 
thousands times, the distribution of the resulted time delay 
was fitted with a Gaussian for inferring the error. 



^ If we assume the excess observed in Group A is real, it can 
be subtracted off in the 30—50 keV flux of bursts in group B by 
scaling the excess derived in group A according to the difference 
in burst peak fluxes averaged over group A and B at 2-10 keV 
band. We therefore subtracted this excess off the 30—50 keV fluxes 
for group B and derived a revised hard X-ray lightcurve in the 
lowest panel of Fig. 4. 

^ We used 1 s-bin lightcurves because for a smaller time interval, 
e.g. 0.5s, the statistics arc insufficient to estimate time lags. 
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where L is the luminosity of the type-I bursts, and C/burst 
and I7disk are the energy densities of the seed photons de- 
rived from the bursts and the disk, respectively. The pro- 
portionality coefficient between them is a. Fig. 7 sketches 
the configuration of the NS-disk-corona system adopted in 
this estimation. The hot corona is located above the in- 
ner accretion disk , which corresponds to western model 
IjWhite et al.lll988l ). where r is the radius and h is the height 
of the corona with respect to the disk. In order to simplify 
the calculatio n, we take the mu lticolor blackbody model 
(T(r)'* oc r"^ (JMitsuda et al.lll984l '). with the innermost tem- 
perature Tin ~ 0.7 keV), and the soft photons dominating 
the innermost disk region. We estimate the disk luminosity 
by approximating it by AnRf^^aT-^. Substituting proper val- 
ues into the above equation (i.e. L = 0.7LEdd, with I/Edd 
representing the Eddington luminosity, and an inner radius 
-Rin = 30 km), [/burst and f/disk/a are shown in Fig. 8 for 
different values of a and the radius r fixed at 25 km. With 
the increase of vertical height of the corona, the reduction 
of (/disk is faster than that of C/burst, and their intersection 
represents the possible height above the disk for which C/disk 
is sufficient for making the system transition into a HSS in 
an outburst. 

To further investigate the possible location of the corona 
while transferring into a HSS, we plot in Fig. 9 the evolution 
of the corona height as a function of the inner disk radius, 
under the condition that f/burst and C/disk are equal. In this 
figure we take an inversely proportional correlation between 
Tlin an d Tin, with a coefficie nt of 14.36 derived from the 
data in lWeng fc Zhang (|201l|). The correlation was derived 
in a form of Tin = 14.36/7?in -I- 0.22. Here a is set to 1, 
and a constant of 0.22 is added to have Tin ~ 0.7 keV at 
-Rin ~ 30 km, consistent with the parameter considerations 
of Fig. 8. Fig. 9 may be regarded as a possible upper limit 
of the vertical height with respect to the underlying disk, 
beyond of which the corona would not be efficiently cooled 
by the disk soft X-rays. Theoretically, a corona could be 
held up stably on top of the disk due to the p ressure balance 
from gravity, magnet and gas. Here we follow iNaravan et al.l 



© 2012 RAS, MNRAS 000, 000-000 



X-ray bursts as a probe of the corona: 4U 1636-536 5 




Figure 7. The NS-disc-corona system adopted in our analysis. 
The blue pancake represents the corona (in order to look at it 
more clearly, we cut a cross section). The colors in the annular 
rings show a cooling accretion disk at larger radii. 



l|l995f) for estimating the possible corona/disk relationship. 
By assuming that a standard viscosity parameter a ~ 0.01, 
and 80 percent energy releasing in the disk propagation to 
the hot flow, a ratio of gas pressure to magnetic pressure ^ 2 
(weak magnetic field) , then a ratio of corona vertical height 
to the radial distance ih/r) can be approximately estimated 
as 0.8 (jNaravan et al.lll995h . Alternatively, a ratio of about 
0. 3 was derived for Sgr A in numerical calculation performed 
bv lManmoto et al.l ([1997]). Therefore, as seen in Fig. 9, a h/r 
value of around 0.3-0.8 would lead to a readily cooled corona 
at a region near the compact object but not once further 
away. We see in Fig. 9 that the disk emission is sufficient for 
cooUng the corona once the latter is located in the vicinity of 
the disk. The corona is more likely cooled by disc emission at 
small inner radii. This is consistent with what is detected in 
outburst of XRBs: changes into a HSS while the disk moves 
towards the compact star. At a larger inner disk radius, the 
evolution shown in Fig. 9 is smoother; in the LHS, the corona 
can be effectively cooled by disk emissions only if located in 
the vicinity of the disk. Those failed outb urst as detected 
in the XRB H 1743-322 l|Chen et aLlbOld ) may have kept 
their corona at least a few km away from the disc during the 
entire outbursts to prevent an effective corona cooling. 



5 DISCUSSION AND SUMMARY 

The shortage of hard X-ray during bursts might be a probe 
of evolution of an XRB outburst. In the LHS of an XRB, the 
hard X-rays can come from a hot corona or a jet. A jet would 
usually show well in the radio band. Quantitatively, the BH 
binaries are ~30 times louder i n radio than atoll s ources 
at similar Eddington luminosity IjMigliari et aLlbOO^ ). The 
correlation between Lj and Lx (the luminosity in radio and 
X-rays, respectively) is quite different between NS and BH 
XRBs: Lj oc Lx in NSs, while Lj oc L^^ in BHs. This dif- 
ference likely indicates that BH s are jet-dominated whil e 
NSs are not jMighari et al.ll2010l: iMigliari fc Fended |2006| ). 
As one example, IMigliari et al.l ( 2010|) concluded that syn- 
chrotron self-Compton (SSC) emission alone, from the post- 
shock jet particles, is not a viable mechanism to explain the 
observed hard X-ray tail of NS 4U 0614+091. This result 




Figure 8. Evolution of C/burst (f^disk/c*) a-s a function of the 
corona height with respect to the disk under different values of 
a, where f/burst arid f7disk are the energy densities of the seed 
photons derived from the bursts and the disk, respectively. The 
intersection shows the possible height of the corona when r is 
fixed at 25 km. 
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Figure 9. The possible height of the corona with respect to the 
disk when C/burst and Cdisk are equal. The vertical axis shows the 
corona height, and the horizontal axis represents the innermost 
radius of the disk during outbursts. 



indicates that hot plasma exists close to the NS in the form 
of a corona, where the hard X-ray emission is associated 
with the accretion disk boundary layers. We thus discuss 
our results mainly in a context of a dominant corona-disk 
radiation. 

To investigate the possible cooling of the corona by soft 
X-ray bursts, we studied 114 bursts embedded in the X- 
ray evolution of 4U 1636-536. We subdivided these bursts 
in two groups, according to the hardness ratio as monitored 
by RXTE/ASM and Swift/BAT. We found that in group B, 
when the source had harder spectrum and therefore higher 
flux at hard X-rays, there is a hint for an anti-correlation 
between the hard and the soft X-rays, similar to what was 
found for IGR 17473-2721. During the bursts, the flux at 
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30-50 keV decreased by roughly 40%, likely indicating that 
the corona was only partially cooled by the ad ditional soft 
X-ray s showers fed by the bursts. As discussed in lChen et alj 
((2013), this may suggest that the corona is not fully covering 
the surface of NS, but that it is perhaps located over the 
disk. 

The hard X-ray shortage lags that of the soft X-rays by 
2.4 ± 1.5 s. Such a time lag is comparable to the one found 
in IGR J17473-2721, and can be regarded as a timescale for 
cooling and recovering of the corona. It is hard to explain 
a corona evolution on a seco nd timescale within a tradi- 
tiona l disk evaporation theorv iJMeyer et al.lll994l : lEsin et al.l 
1 19971 : iLiu et all I2OO7I : iFrank et al.ll2002l ). In these models, 
the corona is thought to be a result of the disk viscosity, and 
therefore its fo rmation has a typi cal timescale of a few days. 
As discussed in lChen et al.l (|2012l ). an alternative model that 
allows such a fast evolutio n is magnetic f ield r e-connection, 
which IZhang et al] (|2000l ). IZhand (|2007l ) and iMaver et al] 
l|2007l) argued for the corona formation in the solar atmo- 
sphere. The magnetic field in the inner accretion disk region 
can be estimated as lO^G IjZhang et al.M200Cll). consistent 
with that inferred from the propeller effect (IWeng fc Zhang 
I2OIII ). However, no characteristic magnetic field -essential 
for re-connection in the disk- is indicated in current mod- 
els, and the disk magnetic field remains unknown. 
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